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ABSTRACT 

We present a study of 581 Hz oscillations observed during a thermonuclear X-ray 
burst from the low mass X-ray binary (LMXB) 4U 1636-54 with the Rossi X-ray Timing 
Explorer (RXTE). This is the first X-ray burst to exhibit both millisecond oscillations 
during the rising phase as well as photospheric radius expansion. We measure an 
oscillation amplitude within 0.1 s of the onset of this burst of 75 ± 17%, that is, almost 
the entire thermal burst flux is modulated near onset. The spectral evolution during the 
rising phase of this burst suggests that the X-ray emitting area on the neutron star was 
increasing, similar to the behavior of bursts from 4U 1728-34 with 363 Hz oscillations 
reported recently. We argue that the combination of large pulsed amplitudes near 
burst onset and the spectral evidence for localized emission during the rise strongly 
supports rotational modulation as the mechanism for the oscillations. We discuss how 
theoretical interpretation of spin modulation amplitudes, pulse profiles and pulse phase 
spectroscopy can provide constraints on the masses and radii of neutron stars. We also 
discuss the implications of these findings for the beat frequency models of kHz X-ray 
variability in LMXB. 
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1. Introduction 

Large amplitude millisecond oscillations have 
now been observed during thermonuclear X-ray 
bursts from six low mass X-ray binary (LMXB) 
systems with the Rossi X-ray Timing Explorer 
(RXTE) (see Strohmayer, Zhang & Swank 1997; 
Smith, Morgan & Bradt 1997, Zhang et al. 1996; 
Swank et al. 1997; and Strohmayer et al. 1997). 
The thermonuclear instability which triggers an 
X-ray burst burns in a few seconds the nuclear 
fuel which has been accumulated on the neutron 
star surface over several hours. This > 10 3 dif- 
ference between the accumulation and burning 
timescales means that it is extremely unlikely 
that the conditions required to trigger the insta- 
bility will be achieved simultaneously over the 
entire stellar surface. This realization, first em- 
phasized by Joss (1978), led to the study of lat- 
eral propagation of the burning instability over 
the neutron star surface (see Fryxell & Woosley 
1982, Nozakura, Ikeuchi k Fujimoto 1984, and 
Bildsten 1995). The subsecond risetimes of ther- 
monuclear X-ray bursts suggests that convection 
plays an important role in the physics of the 
burning front propagation, especially in the low 
accretion rate regime which leads to large igni- 
tion columns (see Bildsten (1998) for a review of 
thermonuclear burning on neutron stars). Bild- 
sten (1995) has shown that pure helium burn- 
ing on neutron star surfaces is in general inho- 
mogeneous, displaying a range of behavior which 
depends on the local accretion rate. Low accre- 
tion rates lead to convectively combustible accre- 
tion columns and standard type I bursts, while 
high accretion rates lead to slower, nonconvective 
propagation which may be manifested in hour 
long flares. These studies emphasize that the 
physics of thermonuclear burning is necessarily 
a multi-dimensional problem and that localized 
burning is to be expected, especially at the onset 
of bursts. 

There is now good evidence that the oscilla- 
tions seen during the rising phase of bursts from 



4U 1728-34 are produced by spin modulation of 
such a localized thermonuclear hotspot on the 
surface of the neutron star, and that the ob- 
served oscillation frequency is a direct measure of 
the neutron star spin frequency (see Strohmayer, 
Zhang & Swank 1997). These observations pro- 
vide the most compelling evidence to date that 
neutron stars in LMXB are rotating with near 
millisecond periods. 

In this Letter we present new burst data from 
the LMXB 4U 1636-54 which provides further ev- 
idence in support of the spin modulation hypoth- 
esis for the millisecond burst oscillations. We 
present data from a thermonuclear burst from 
4U 1636-54 which reveals a strong transient os- 
cillation during the burst rise at 1.723 ms with an 
initial amplitude of 75 ±17%. We also discuss the 
implications of these findings for the spin mod- 
ulation interpretation and how they can be used 
to place constraints on the mass and radius of 
neutron stars. Finally, we discuss some implica- 
tions of our observations for the current theories 
of kilohertz quasiperiodic oscillations (QPO) in 
LMXB. 

2. Observations and data description 

Zhang et al. (1996) reported the discovery of 
581 Hz oscillations during several thermonuclear 
X-ray bursts from 4U 1636-54. Here we focus 
on a single burst detected from this source with 
RXTE at 22:39:22 UTC on December 28, 1996. 
This burst is unique in that it reveals strong os- 
cillations at 581 Hz prior to the onset of photo- 
spheric radius expansion. Oscillations are not de- 
tected during the photospheric expansion phase, 
but reappear with much weaker amplitude after 
photospheric touchdown in a manner very sim- 
ilar to that reported for bursts from 4U1728-34 
and KS 1731-26 (Strohmayer et al. 1997; and 
Smith, Morgan, & Bradt 1997). 

For this burst we obtained event mode data 
with 125 /is (1/8192 s) time resolution and 64 
spectral channels across the 2-100 keV PC A 
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bandpass. This data mode employs a pair of data 
buffers which are alternately read out after a se- 
lectable accumulation interval. Each buffer can 
store ps 16,000 events. In the case that the source 
countrate is high enough to fill the buffer before 
the end of an accumulation interval gaps in the 
data can occur. For the event data shown here 
we did have short gaps when the burst flux was 
high, but telemetry constraints preclude selecting 
a shorter accumulation interval. For these bursts 
we also had 125 ^s, 2-90 keV binnned timing data 
which does not suffer from gaps. 

3. Oscillations during the rise of a radius 
expansion burst 

At the time of the burst in question the persis- 
tent countrate from 4U1636-54 in the total PCA 
bandpass was about 2100 cts/sec. 4U1636-54 is 
classified as an atoll source and it was in the 
so called banana branch during these observa- 
tions (see Hasinger & van der Klis 1989). We 
will present more detailed analysis of the accre- 
tion driven flux and properties of all the X-ray 
bursts observed in a subsequent paper. Figure 
1 shows the 2 - 90 keV PCA lightcurve (top) 
and the hardness ratio (6 - 18 keV) / (2 - 6 keV) 
at 1/16 s resolution for this burst. The drop in 
the hardness ratio at about 1.5 s is a clear in- 
dication of photospheric radius expansion. This 
burst had a peak flux in the 2-60 keV band of 
7.2 x 10~ 8 ergs cm -2 s _1 and a fluence in the 
same band of 3.9 x 10 ergs cm -2 . At the sys- 
tem's probable distance of 6.5 kpc (see Lawrence 
et al. 1983) these correspond to 3.6 x 10 38 ergs 
s~ l and 2.0 x 10 39 ergs. We used the 2 - 90 keV, 
125 /us binned data to compute a dynamic power 
spectrum by calculating power spectra from 2 s 
of data every 1/8 s through the burst. Since 
the data overlap in this analysis the individual 
power spectra are not independent, however, this 
procedure provides a straightforward way to in- 
vestigate the time evolution of the oscillations. 
The resulting dynamic power spectrum is pre- 
sented in figure 2, which shows contours of con- 



stant power spectral amplitude along with the 
burst profile superposed. This analysis reveals 
the presence of strong oscillations at about 580 
Hz on the rising edge of the burst as well as os- 
cillations after photospheric touchdown at a fre- 
quency about 1 Hz higher. As we discuss below, 
the oscillations during the rise are transient, last- 
ing about 0.25 s, this introduces a characteristic 
width to the peak in the power spectrum which 
results in the broader distribution in frequency of 
the power spectral contours near the burst rise. 
We emphasize that the contours centered near 
580.5 Hz at about 2 s are not indicative of a real 
drop in the oscillation frequency at the burst on- 
set. However, the change in frequency from burst 
rise to post-peak is real and is not caused by fi- 
nite sampling time effects or the transient nature 
of the oscillations. This trend for the oscilla- 
tion frequency to increase with time during some 
bursts has now been reported in all sources with 
millisecond burst oscillations except KS1731-26, 
from which only a single burst with oscillations 
has so far been detected (see Strohmayer et al. 
1997; Smith, Morgan, & Bradt 1997). 

To further characterize the oscillations seen 
during the rising phase we performed an epoch 
folding period search on the rising interval. Fig- 
ure 3 shows the lightcurve of the burst, focusing 
on the rising portion. The interval used for this 
analysis is denoted by the dashed lines in figure 
3. The epoch folding reveals a strong oscillation 
at 1.723 ms. We folded the data on this period to 
determine the pulse profile and amplitude. The 
resulting pulse profile is shown in figure 4a. We 
fit a model with a constant countrate plus a si- 
nusoid, A + Bsm{2^t/P + 4>), to the pulse profile 
and find an amplitude B/A = 30 ± 4% with an 
acceptable fit. Since there is no evidence that 
the persistent flux component is oscillating at 
the frequency observed during the burst, to de- 
termine the pulsed amplitude we first subtracted 
from the folded lightcurve the average persistent 
countrate calculated from a 20 s interval prior to 
the burst. We note that there are no detectable 
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oscillations during the interval when the drop in 
hardness ratio at about 1.5 s (see figure 1) in- 
dicates the onset of photospheric radius expan- 
sion. This is consistent with the idea that the 
expansion of the photosphere destroys the co- 
herent oscillations, probably both by scattering 
away the amplitude as well as changing the fre- 
quency quickly via the expansion of the surface 
layers. After photospheric touchdown (see figure 
1 at about 3.5 s) the oscillations are again de- 
tected, but with slightly higher frequency and a 
much lower amplitude of about 5.7 % in the 2 - 
60 keV bandpass. 

The presence of large amplitude oscillations 
very near the onset of this burst is similar to the 
behavior of bursts from 4U 1728-34 reported by 
Strohmayer, Zhang & Swank (1997). We found 
even higher oscillation amplitudes by folding a 
62 ms interval close to the onset of the burst. 
The folded interval is denoted by the solid verti- 
cal lines in figure 3. The background subtracted 
folded lightcurve obtained from this 62 ms in- 
terval is shown in figure 4b. We again fit the 
sinusoidal model to this lightcurve and found an 
amplitude of 75 ± 17% for this interval. This is 
easily the largest amplitude of the 581 Hz os- 
cillation seen during this burst. Such a large 
amplitude-almost the entire thermal burst flux 
is modulated-only about 0.1 s after the onset of 
the burst is strong evidence that we are seeing a 
rotationally modulated, thermonuclear hotspot. 

Strohmayer, Zhang & Swank (1997) reported 
spectral evidence for an increasing X-ray emit- 
ting area during burst rise using bursts from 
4U 1728-34 observed with RXTE. We performed 
similar spectral analysis in order to determine 
whether the same process is plausible for this 
burst as an explanation for the strong oscillations 
seen during the rise. Using the 125 /xs, 64 channel 
data we accumulated spectra through the burst 
on either 0.0625, 0.125, 0.25 or 0.5 s intervals 
depending on the countrate. We also accumu- 
lated a 20 s interval just prior to the burst for 
use as a background estimate. We fit the back- 



ground subtracted spectra to a blackbody model 
to determine temperatures and fluxes as a func- 
tion of time through the burst. We created a 
plot of F^ 4 /kT BB versus F bo i, where F boi is the 
bolometric flux and ¥Tbb is the blackbody tem- 
perature in keV. For a true blackbody the quan- 
tity F^/kTsB oc A l / A /d 1 / 2 , where A is the X- 
ray emitting area and d is the distance to the 
source. We found results very similar to those 
reported by Strohmayer, Zhang & Swank (1997) 
for bursts from 4U 1728-34, that is, F^f /kT BB 
increases during the rise. This is strong evidence 
for an expanding X-ray emission region and thus 
supports the spin modulation hypothesis in this 
burst from 4U 1636-54. We emphasize that this 
result suggests an increase in the X-ray emitting 
area prior to the onset of photospheric radius ex- 
pansion. The onset of photospheric radius expan- 
sion is reflected in the drop in hardness ratio at 
about 1.5 s in figure 1. It also corresponds to 
an increase in the emitting surface area, but it 
is caused by radial expansion of the photosphere 
due to radiation forces and not the lateral propa- 
gation of the burning instability. The blackbody 
spectral analysis indicates that the inferred emit- 
ting areas were approximately the same both im- 
mediately before and after the photospheric ra- 
dius expansion. This result combined with the 
large difference in oscillation amplitudes before 
and after photospheric expansion suggests that 
the thermonuclear flash covered the entire sur- 
face of the neutron star. 

4. Discussion 

If the millisecond oscillations seen during bursts 
are in fact due to spin modulation, then de- 
tailed study and modelling of the oscillation am- 
plitudes, pulse profiles and spectral variability 
with pulse phase during X-ray bursts can provide 
a wealth of information on the mass and radius 
of the neutron star. For example, the maximum 
modulation amplitude that can be obtained from 
a hotspot of a given angular size on a rotating 
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neutron star is set by the strength of general rel- 
ativistic light bending. For the case that rotation 
is not rapid enough to substantially distort the 
exterior spacetime from the Schwarzchild space- 
time, and this is the case even for spin periods of 
a few milliseconds (Lamb & Miller 1995), then 
the maximum amplitude depends only on the 
compactness of the neutron star, that is, the ratio 
of stellar mass to radius GM/c 2 R. Stars which 
are more compact produce lower amplitudes due 
to flattening of the pulse by light bending (see 
Pechenick, Ftaclas & Cohen 1983; Strohmayer 
1992; and Miller k Lamb 1997). Since the intrin- 
sic rotational modulation amplitude can only be 
decreased by other effects such as photon scat- 
tering (see Miller, Lamb & Psaltis 1997; Brain- 
erd & Lamb 1987; and Kylafis & Phinney 1989) 
or the viewing geometry of the spot, the maxi- 
mum observed oscillation amplitude represents a 
lower limit to the intrinsic amplitude. Thus an 
observed amplitude can be used to place an up- 
per limit on the compactness of the neutron star, 
that is, if the star were more compact than some 
limit it would not be able to produce a modula- 
tion amplitude as large as that observed. 

In principle, stellar rotation will also play a 
role in the observed properties of spin modula- 
tion pulsations. For example, assuming the os- 
cillation frequency of 581 Hz represents the spin 
frequency of the neutron star in 4U 1636-54, then 
for a 10 km radius neutron star the spin velocity 
is v sp i n /c = 2irv sp i n R ss 0.12 at the rotational 
equator. The motion of the hotspot produces a 
Doppler shift of magnitude AE/E ~ v sp i n /c = 
0.12, thus the observed spectrum is a function of 
pulse phase (see Chen & Shaham 1989). Mea- 
surement of a pulse phase dependent Doppler 
shift in the X-ray spectrum would provide addi- 
tional evidence supporting the spin modulation 
model and would also provide a means of con- 
straining the neutron star radius. The rotation- 
ally induced velocity also produces an aberra- 
tion which results in asymmetric pulses, thus the 
pulse shapes also contain information on the spin 



velocity and therefore the stellar radius (Chen & 
Shaham 1989). The component of the spin veloc- 
ity along the line of site is proportional to cos 9, 
where 6 is the lattitude of the hotspot measured 
with respect to the rotational equator. The mod- 
ulation amplitude also depends on the lattitude 
of the hotspot, as spots near the rotational poles 
produce smaller amplitudes than those at the 
equator. Thus we expect a correlation between 
the observed oscillation amplitude and the size of 
any pulse phase dependent Doppler shift. Dec- 
tection of such a correlation in a sample of bursts 
would definitively confirm the rotational modula- 
tion model in our opinion. We will present calcu- 
lations of mass - radius constraints for 4U 1636- 
54 and 4U 1728-34 including the effects of light 
bending, rotation and angle dependent emission, 
based on the observed properties of burst oscil- 
lations as well as spectroscopy of Eddington lim- 
ited bursts in a subsequent paper. 

Several models for the kilohertz QPO seen in 
13 LMXB systems (see van der Klis 1997 for a re- 
cent review) invoke some sort of beat-frequency 
interpretation for the twin kHz peaks seen in 
many of the sources (see Miller, Lamb & Psaltis 
1997; Strohmayer et al. 1996). So far only 
in 4U 1728-34 does the frequency difference be- 
tween the twin kHz peaks match the frequency 
observed during X-ray bursts (see Strohmayer et 
al. 1996). In two other sources (KS 1731-26 and 
4U 1636-54) the separation of the twin kHz peaks 
is closer to 1/2 the frequency of oscillations ob- 
served during bursts (see Wijnands & van der 
Klis 1997; and Zhang et al. 1996). For exam- 
ple, Wijnands et al. (1997) report a frequency 
difference for the twin kHz QPO in 4U1636-54 of 
276±10 Hz. The effort to reconcile these observa- 
tions with a beat-frequency interpretation has led 
to speculation that the oscillation frequency ob- 
served during bursts may sometimes be twice the 
spin frequency of the star, although in 4U1636- 
54 the difference frequency appears to be a bit 
less than 1/2 the burst oscillation frequency, and 
in some Z sources the frequency difference is not 
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constant (van der Klis et al. 1997). If this sce- 
nario is correct it implies the existence of two 
antipodal spots on the neutron star surface dur- 
ing X-ray bursts. In addition to the daunting 
requirement of initiating the thermonuclear flash 
nearly simultaneously on opposite sides of the 
neutron star, the observation of large oscillation 
amplitudes shortly after burst onset in 4U1636- 
54 places severe constraints on the two hotspot 
scenario. To see this one can ask the following 
question. What maximum amplitude can be pro- 
duced by a star with antipodal hotspots? For 
two antipodal spots light bending strongly con- 
strains the amplitudes that can be achieved (see 
Pechenick, Ftaclas & Cohen 1983). Calculations 
using the Schwarzchild spacetime and isotropic 
emission from the stellar surface indicate that 
even a neutron star with an implausibly small 
compactness of M/R = 0.1, recall that rotational 
modulation amplitude increases with decreasing 
compactness, can only achieve a maximum am- 
plitude of about 30 %, whereas we measured an 
amplitude of 75% from the burst described here. 
We note that with M/R = 0.1 a 1.4 M sun neu- 
tron star would have a radius of 21 km. This is 
far stiffer than any neutron star equation of state 
that we are familiar with. Rapid rotation could 
in principle modify this result, but at the mod- 
est inferred spin period of 290.5 Hz for the two 
spots to produce the observed 581 Hz frequency 
one would still require implausibly large neutron 
star radii to make a significant rotational correc- 
tion to the amplitude. Another process which 
can increase the amplitude is beaming of radia- 
tion at the stellar surface. For example, to the 
extent that electron scattering is the dominant 
opacity process in neutron star atmospheres then 
one should expect a specific intensity distribu- 
tion which approximates that from a grey atmo- 
sphere. Such a distribution is proportional to 
cos 5 + 2/3, where 5 is the angle from the nor- 
mal to the stellar surface (see Mihalas 1978), so 
this will modestly increase the rotational mod- 
ulation amplitude. In addition, it is likely that 



the emergent spectrum will also be a weak func- 
tion of 5 (Miller & Lamb 1997). With more de- 
tailed modelling of the emission from hotspots of 
finite size it will be possible to place strong con- 
straints on the antipodal hotspot interpretation 
for the burst oscillations. Finally, if further anal- 
ysis such as pulse phase spectroscopy continues 
to support the interpretation of the burst oscilla- 
tion frequency as the neutron star spin frequency 
in 4U1636-54, then kHz QPO frequency separa- 
tions near 1/2 the spin frequency would suggest 
that the beat frequency interpretation may be 
untenable. 
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5. Figure Captions 

Figure 1: 2-90 keV lightcurve and hardness ratio 
plot for the burst. The hardness ratio is defined 
as the ratio of counts in the 6-18 keV range to 
that in the 2-6 keV range. The drop in hardness 
ratio at about 1.5 s indicates the beginning of 
the radius expansion phase. The dashed vertical 
lines denote the time interval used in the epoch 
folding analysis. 

Figure 2: Representation of the dynamic power 
spectrum for the burst shown in figure 1. The 
contours are levels of constant power spectral am- 
plitude. The power spectra were calculated from 
2 s intervals with the center of each succesive in- 
terval shifted by 1/8 s. The countrate profile of 
burst 1 is also shown. The burst profile is aligned 
to the center of each data interval used to calcu- 
late the power spectra. 

Figure 3: Lightcurve of the burst focusing on the 
rise. The interval denoted by dashed lines is that 
used for the folding analysis shown in figure 4a, 
while the solid vertical lines denote the 62.5 ms 
interval used for the folding analysis displayed in 
figure 4b. 

Figure 4: (a) Folded, background subtracted 
lightcurve using the 0.25 second interval denoted 
with dashed vertcal lines in figure 3. The best 
fitting sinusoidal model is also shown. The aver- 
age amplitude during this interval was 30 ± 4%. 
(b) The same results are shown for the 62.5 ms 
interval shown in figure 3. For this interval the 
amplitude was 75 ± 17%. 
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